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Background: Evidence increasingly supports that prostate cancer is initiated by the malignant transformation of stem cells (SCs). 
Furthermore, many SC-signalling pathways are shown to be shared in prostate cancer. Therefore, we planned transcriptome 
characterisation of adult prostate SCs as a strategy to consider new targets for cancer treatment. 

Methods: Intuitive pathway analysis was used for putative target discovery in 12 matched selections of human prostate SCs, 
transiently amplifying cells and terminally differentiated cells. These were pooled into three groups according to the stage 
of differentiation for mRNA microarray analysis. Targets identified were validated using uncultured primary tissue (n = 12), 
functional models of prostate cancer and a tissue microarray consisting of benign (n = 42) and malignant prostate (n = 223). 

Results: A deficiency in class 1 UDP glucuronosyltransferase (UGT) enzymes (UGT1A) was identified in prostate SCs, which are 
involved in androgen catabolism. Class 1 UGT enzyme expression was also downregulated in cancer SCs and during progression 
to metastatic castration-resistant prostate cancer (CRPC). Reduction of UGT1 A expression in vitro was seen to improve cell survival 
and increase androgen receptor (AR) activity, as shown by upregulation of prostate-specific antigen expression. 

Interpretation: Inactivation of intracellular androgen catabolism represents a novel mechanism to maintain AR activity 
during CRPC. 



Mechanistic studies of tumorigenesis point towards stem cells 
(SCs) being the cell of origin in prostate cancer initiation, 
suggesting stem pathways have a crucial role in the early stages 
of cancer development (Gu et al, 2007; Goldstein et al, 2010; 
Lawson et al, 2010). Furthermore, known stem pathways have 
been seen to be aberrantly expressed as the prostate cancer 
progresses (Linn et al, 2010; Jia et al, 2011; Kregel et al, 2013). 
Characterisation of the genes at play within prostate SCs, and 
characterising these same markers in the malignant setting may 
identif)^ unique therapeutic opportunities for the treatment of 
prostate cancer. 



We have recently demonstrated the presence of the androgen 
receptor (AR) within normal human prostate SCs (Williamson 
et al, 2012), so identification of specific novel AR pathway 
processes within these cells will be of particular importance. Such 
mechanisms may explain the maintenance of androgen 'addiction' 
in the malignant setting during cancer progression despite 
androgen deprivation, as it is established that the AR pathway 
has diverse roles depending on the degree of prostate differentia- 
tion (Gregory et al, 1998; Buchanan et al, 2001; Mousses et al, 
2001; Chen et al, 2004; Taplin and Balk, 2004; Attard et al, 2009; 
Bonkhoff and Berges, 2010; Lee et al, 2012). 
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Human prostate SCs can be enriched by their expressions of 
a2i^i-integrins and CD 133 (Richardson et al, 2004). In this work, 
we describe gene expression signatures associated with SCs 
((X2p^^ CD133 + ^^), transiently ampUfying (TAP) cells (asj^f 
CD 133"^^) and terminally differentiated (TD) cells {(X2p\^'^ 
CD133~^^) following the intuitive pathway analysis of transcrip- 
tome characterisation of human prostates. Within normal prostate 
SCs, we discovered a deficiency in class 1 UDP (uridine 
5'-diphospho)-glucuronosyltransferase; UGT) enzymes (UGTAl), 
which are regulators of androgen catabolism. Our functional 
in vitro models of prostate cancer demonstrated that UGTIA 
knock down leads to upregulation of prostate-specific antigen 
(PSA) along with increased cell survival. Furthermore, reduced 
UGTIA expression was shown to be a feature of clinical prostate 
cancer that was associated with a poorer prognosis. 



MATERIALS AND METHODS 



Tissue collection and gene expression array analysis. Human 
prostate samples were collected from a total of 50 patients and 
processed, as previously described (WilUamson et at, 2012), to 
yield asj^™ CD133 + ^^ cells (stem progenitors), a2i^™CD133 "^^ 
(transiently amplifying progenitors) and (X2P\^^ CD133~^^ 
(differentiated) epithelial cells without culturing (Richardson 
et al 2004; Heer et a/, 2006; Williamson et al 2012). These 
immediately underwent RNA isolation (micro RNeasy, Qiagen, 
Manchester, UK), with only the highest quality samples being 
taken forward for analysis. To ensure the highest quality of data, 
before being taken forward for study, samples must have 
demonstrated an RNA integrity number >8.0 in each selected 
fraction. This approach yielded 12 high-quality matched patient 
samples to be taken forward for further analysis (Supplementary 
Table 1). These were pooled using a validated methodology 
to account for biological variance (Kendziorski et al, 2005; 
Zhang et al, 2007) creating three pooled samples of a2i^™ CD133^^, 
a2i^f CD133"^ and a2p\^^ CD133"^. These samples were then 
investigated on Afiymetrix (High Wycombe, UK) U133 plus 2.0 chips. 
Detailed methods and analyses used are outlined in depth in the 
Supplementary Information. 

Quantitative real-time PGR analysis. MessageBOOSTER cDNA 
synthesis amplification kit (Epicentre Biotechnologies, Madison, 
WI, USA) was employed for prospective validation of gene 
signatures identified in the microarray on the unpooled 
individual sorted samples. Real-time PGR (Applied Biosystems 
(Life Technologies, Paisley, UK) 7900HT) was performed on 
subsequent cDNA libraries using SYBR Green (Invitrogen, Life 
Technologies, Paisley, UK) and cDNA-specific primers sets 
detailed in Supplementary Table 2. All data were normalised to 
the housekeeping gene GAPDH. 

Maintenance of prostate cancer cell lines. Human prostate 
benign cell line BPH-1 in addition to human prostate cancer cell 
lines LNCaP and VCaP (American Type Culture Collection 
(ATCC), Manassas, VA, USA) were maintained in RPMI 1640 
medium (Sigma, Gillingham, UK) containing 10% fetal calf serum 
(Sigma) and 2mM L-glutamine (Sigma), referred to as complete 
medium. The isogenic androgen -independent cell line LNCaP-AI, 
used to mimic castration resistant prostate cancer (CRPC; 
Halkidou et al, 2003), was maintained in steroid-depleted media; 
RPMI 1640 medium (Sigma) containing 10% heat- inactivated 
dextran charcoal-treated fetal calf serum (Hyclone, Fisher 
Scientific, Loughborough, UK) and 2 mM L-glutamine (Sigma). 

Small-interfering RNA knock down. Cells were seeded in 
six- well plates before being transfected with either non- silencing 
small-interfering RNA (siRNA) or equivalent concentration 



UGTIA siRNA (sequences are listed in Supplementary 
Information) using Lipofectamine RNAiMAX (Invitrogen). Addi- 
tional methods, including western blot, are available in 
Supplementary Information. 

Clonogenicity assays. Colony-forming assays were performed by 
plating 2000 viable cells in individual wells of a six- well dish for 14 
days, fixing in Carnoy s frxative (Sigma) and counterstaining with 1% 
Crystal Violet (Sigma). Colony counting was performed by ColCount 
automated colony counter (Oxford Optronix, Milton, UK). 

Tissue microarray staining and scoring. Class 1 UGT enzyme 
staining was carried out on a tissue microarray consisting of benign 
(n = 42) and malignant prostate (n = 223) material derived from 
radical prostatectomy or trans-urethral resection of the prostate 
(Heer et al, 2004). Sections were deparaffinised and rehydrated 
before undergoing antigen retrieval under pressure in 0.01m 
sodium citrate buffer and were treated for endogenous peroxidase 
activity. Sections were stained with 1 : 100 UGTIA antibody 
(Sc-25847, Santa Cruz Biotechnology, Heidelberg, Germany) and 
the Menapath two-step HRP universal polymer kit (A. Menarini 
Diagnostics Ltd., Winnersh- Wokingham, UK), according to the 
manufacturer's recommendations. Immunostaining was classified as 
absent (0), weak (1), moderate (2) or strong (3) by a consultant 
uropathologist (Mantilla). Staining was correlated to clinical 
parameters using Kruskal-Wallis test, Mann-Whitney L^-test, 
Log- rank analysis and Cox regression analysis. Where appropriate, 
P<0.05 was considered statistically significant. 



RESULTS 



Post-hybridisation quality assurance of microarray data series. 

Quality assessment of samples confirmed 3'-end/-middle (37M) 
ratios for housekeeping genes below the required threshold of 3, 
which confirms linear amplification of samples before running of 
samples and was in keeping with the high-quality threshold set for 
the samples utilised (Supplementary Figure S3). Investigation of 
the top genes up- and downregulated between the epithelial 
fractions confirmed CD 133 as the top ranked differentially 
expressed transcript in the comparison with the a2j^™ CD133^^ 
and (X2p^^ CD133"^^ fractions (17.88-fold increase in the asj^™ 
CD133^^^ fraction). Distinct gene signatures were observed for 
each of the transcript profiles for prostate SCs, TAPs and TD cells, 
and prospective validation of these transcript patterns during 
differentiation was confirmed using CD 133 and p63 real-time PGR 
(Supplementary Figure S3), with patterns of expression for key 
markers of prostate differentiation (AR, KLK2, KLK3 and 
TMPRSS2) matching those shown previously (Williamson et al, 
2012). 

Differential androgen metabolism and signalling is observed in 
SCs, CSCs and throughout disease progression. Pathway 
ontology analysis was undertaken to characterise groups of genes 
altered during prostate differentiation based on their interacting 
biological function. Pathway analyses comparing the a2Pi 
CD 133 + ^^ fraction of SCs with the (X2p^^ CD 133"^^ fraction of 
early differentiation-committed cells (TAPs) highlighted 'testoster- 
one biosynthesis, metabolism' as one of the top ranked 
differentially expressed pathways (P< 0.001; Figure 1). Detailed 
interrogation of the androgen pathways revealed no change in gene 
expression within the anabolic component of testosterone 
biosynthesis 

during normal prostate differentiation (Supplementary Figure S4). 
However, the catabolic pathways showed a significant difference 
but only in the UGTlA-enzyme family. Class 1 UGT enzyme 
expression showed 4- to 10-fold downregulation in normal SCs 
compared with more differentiated cells, depending upon the 
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A Top 1 0 statistically significant pathways identified by metacore 

Cell adhesion and ECM remodelling ^^B^M 
Testosterone biosynthesis and metabolism | | 
Development of TGF-j9-induced EMT via SMADs | 
Cell adhesion plasmin signalling | 
Development of WNT-signalling pathway I 
G-protein signalling 1 
Cell adhesion gap junctions 1 
Signal transduction calcium signalling | 
Oestradiol metabolism | 
TGF-p regulation of cell proliferation | 
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Figure 1 . Pathway analysis of prostate hierarchy. (A) (GeneGo Inc., St Joseph, Ml, USA) analysis comparing genes two-fold different between the 
oczP^^ CD133^^^ and oczP^^ CD133~^^ fractions. Top 10 pathways are ranked according to their statistical significance (inverse log P-value). (B) 
Metacore-annotated pathway map (GeneGo Inc.) for testosterone biosynthesis and metabolism comparing the oczP^^ CD133^^^ with the a2iSV' 
CD133~^^ fraction. Blue thermometers indicate downregulation of genes in the oczP^^ CD133^^^ fraction. 



different probe sets for this family in the array (Figure IB and 
Supplementary Figure S5). UDP glucuronosyltransferase enzymes 
remove hydrophobic molecules, such as testosterone, from cells by 
glucuronidation. Previous studies on testosterone glucuronida- 
tion have concentrated on the class 2 enzymes in the prostate and 
in prostate cancer, specifically UGT2B7, UGT2B15 and 
UGT2B17 (Belanger et al, 2003; Chouinard et al, 2007); however, 
these were seen to be unaltered between the three epithelial 
fractions of this study (Supplementary Figure S5). Of note, the 
specific probe set for UGT1A8 was the third most downregulated 
probe set (10-fold, P< 0.001) in the comparison between the 
a2pf^ CD133 + ^^ and the asj^™ CD133"^^ fraction. Assessment 
of individual UGTIA enzymes using isoform-specific real-time 
PGR primer sets in unpooled individual human prostate samples 
confirmed UGT1A8 as the highest expressed enzyme in the 
(X2p\^'^ GD133"^^ TD fraction, and that none of the UGTAl 
enzymes were detectable in the asj^™ CD133 + ^^ SG fraction 
(Figure 2A). Meta-analyses of UGTAl enzyme expression in 
studies of prostate SGs and cancer SGs (GSGs; Birnie et at, 2008) 
corroborated these findings, showing decreased expression in 
both SGs and GSGs compared with differentiated cells (P<0.01, 
gene set enrichment analysis (GSEA); Figure 2B). Moreover, 
assessment of UGTAl enzymes in clinical transcriptome data 
from cases of prostate cancer progression (Varambally et al, 
2005) confirmed further reduced expression in GRPG compared 
with early prostate cancer (P<0.01, GSEA; Figure 2G). Reduced 



UGTIA expression in progenitor cells suggests a potential 
mechanism for the accumulation of androgens within prostate 
cancer cells during androgen deprivation therapy, and UGTIA 
expression remains depressed in GRPG. 

Reduction in UGTIA expression is seen in cell line models of 
CRPC and leads to increased AR activity. The functional effects 
of reduced UGTAl enzymes expression were explored using 
functional models of prostate cancer. Glass 1 UGT protein 
expression was confirmed in the benign prostatic cell line BPH-1 
and the androgen- dependent cell line LNGaP; however, reduced 
levels of expression for UGTIA were seen in both the isogenic 
LNGaP- AI and castration-resistant derived line VGaP (Figure 3A). 
To explore whether reduced UGTIA expression was a functional 
feature of prostate cancer progression, LNGaP cells were 
transfected with siRNA targeting the UGTAl enzymes in standard 
steroid-containing medium. Following UGTIA knock down, 
increased AR activity was evident by increased PSA expression 
despite no change in AR expression compared with scrambled 
siRNA control (Figure 3B). Furthermore, following UGTIA knock 
down LNGaP cells were seen to have a two -fold increase in 
clonogenicity compared with scrambled siRNA controls (P<0.05; 
Figure 3G). These data confirm that AR is more active in cells, 
which is probably because of the presence of relatively higher 
intracellular levels of steroids expected with lower levels of 
UGTIA. 
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Figure 2. Expression of UGTA1 enzymes in the prostate 
differentiation and prostate cancer progression. Error bars represent 
s.e. of the mean. (A) Expression of individual UGTA1 enzymes following 
real-time PGR normalised to GAPDH (n = 6). (B) The GSEA for UGT1A 
comparisons in stem cell (SG) and cancer SG (GSG) vs terminally 
differentiated (TD) normal and terminally differentiated cancer (GTD) 
cells, respectively. Barcode plots are graphical representations of the 
Wilcox gene-set test results using ranks. The ranks of the sorted 
log-fold changes between conditions run along the x axis from high 
(pink) to low (green). Black lines represent matches to the UGTA1 
enzymes. P-values are given in the tables below for the alternative 
hypothesis that the UGTA1 enzymes tend to be (i) upregulated in the 
data set or (ii) down regulated. (C) The GSEA for comparisons of UGT1 A 
expression in benign prostate (Ben.), primary prostate cancer (1°) and 
castrate-resistant metastatic cancer (Mets). 
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Figure 3. Functional characterisation of altered UGT1 A expression in 
prostate cancer. (A) Representative western blot analysis measuring 
expression of UGT1A. (B) Representative western blot analysis of 
UGT1 A, AR and PSA following knock down with UGT1 A siRNA. (C) A 
2D colony-forming efficiency of LNGaP in complete medium following 
knock down with UGT1A siRNA (n = 3; P<0.05). 



Class 1 UGT enzyme expression is downregulated in clinical 
prostate cancer progression. Decreasing levels of UGTIA 
immunostaining was associated with increasingly aggressive 
disease, as determined by D'Amico prostate cancer risk 
(based on Gleason grade, stage and presenting PSA; P = 0.04, 
Kruskal-Wallis test; Figure 4A and B). Low levels of UGTIA 
expression were also associated with increased serum PSA levels 
upon clinical presentation (P = 0.003, Mann- Whitney U-test; 
Figure 4C), consistent with reduced UGTIA expression increasing 
intracellular testosterone bioavailability and AR activity. Interest- 
ingly, despite patients with low UGTIA expression being 
associated with higher-risk disease, they initially responded much 
better to hormonal therapy than those with high UGTIA 
expression (P = 0.03, Mann-Whitney (7-test; Figure 4D), suggest- 
ing a more pronounced ligand- dependant AR activity in tumours 
with low UGTIA expression. However, this initial responsiveness 
to hormones is transient as reduced UGTIA expression correlated 
with a significant reduction in patient survival (P = 0.025, 
Log-rank analysis; Figure 4E). In particular, this survival effect 



was pronounced in patients with metastatic disease (P = 0.038, 
Log- rank analysis; Figure 4F). Interestingly, time to CRPC was not 
significantly different between high- and low-UGT groups, and it 
may be that different hormone-escape pathways are established in 
these two groups during progression to death from the advent of 
CRPC. Although Cox multivariate analysis showed that UGTIA 
expression did not confer additional prognostic indication of 
survival independent of grade and stage, we nevertheless show that 
reduced UGTIA expression correlates with an increased AR 
signalling, supporting its role as a mechanism of castration 
resistance. 



DISCUSSION 



There is growing evidence that it is the SCs in the ageing prostate 
that acquire the required mutations to become the cell of origin for 
malignant transformation and that prostate tumours are main- 
tained by a hierarchal CSC model (Collins et at, 2005; Gu et at. 
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Figure 4. Expression of UGT1 A within a tissue microarray of prostate cancer (n = 265; benign n = 42, nnalignant n = 223). Error bars represent s.e. 
of the mean. (A) Representative images of UGT1 A staining in prostate tissue. (B) Staining intensity for UGT1 A for samples stratified for D'Amico 
risk. (C) Presenting PSA levels according to staining intensity for UGT1 A. (D) Prostate-specific antigen nadir levels following hormonal therapy 
according to UGT1 A expression (E) A Kaplan-Meier plot demonstrating survival according to UGT1 A staining intensities of 0-1 (blue line) and 2-3 
(green line). (F) Kaplan-Meier plots of prostate cancer survival in patients without metastases (left hand plot) or with metastases (right hand plot) 
with UGT1A staining intensity of 0-1 (blue line) or 2-3 (green line). 



2007; Goldstein et al 2010; Lawson et al 2010; Williamson et al 
2012). However, the relevance of the signalling pathways 
maintaining these cells types during progression of the disease 
into lethal CRPC remains poorly understood. A transcriptome- 
based target discovery approach, looking at the pathways within 
the SC hierarchy and their role in the malignant setting would be 
of interest. This approach was explored in this paper and revealed 
novel target pathways. 

There is an ongoing debate over the nature of the SC 
organisation of the prostate epithelium. Most models of prostate 
homoeostasis have advocated the presence of undifferentiated SCs 
within the basal compartment of the prostate (English et al, 1987; 



Isaacs and Coffey, 1989; Richardson et a/, 2004). In particular, a 
single basally derived SC has been shown to regenerate a fully 
differentiated prostate epitheUum (Leong et at, 2008). The location 
of the SC niche, however, has been questioned following the 
identification of luminally derived castration-resistant Nkx3.1 cells 
(CARNs) (Germann et at, 2012), which are also capable of 
reconstituting the prostate epithelium from a single cell. An 
alternative model comes from mouse lineage-tracing experiments 
that show basally derived multipotent SCs give rise to unipotent 
basal and luminal progenitors during postnatal development; 
however, in the adult gland, the epithelium is maintained 
by unipotent progenitors specific for basal and luminal layers 
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(Choi et al, 2012; Ousset et al, 2012). Although mice models 
remain an important tool for the prostate modelling of 
tumorigenesis and differentiation, studies of human material 
remains the gold standard and continue to support a basal SC in 
the human prostate. Recent lineage tracking in human tissue 
demonstrated that all epithelial cells share a clonal origin 
(Blackwood et al, 2011) and most published work supports a 
human basal SC (Richardson et al, 2004; Goldstein et al, 2008). 

One of the issues faced when investigating small numbers of 
samples (<12) with microarray technology is the potential for 
outlier patient- specific differences masking more critical under- 
lying population differences between fractions (Kendziorski et al, 
2005; Zhang et al, 2007). To minimise the potential for this in our 
study, we opted to investigate more the transcriptome of the 
prostate as a target discovery platform using a validated pooling 
strategy. This pooling strategy has been demonstrated to account 
for biological variance within populations when ^12 cases are 
included and has been validated in gene array studies (Kendziorski 
et al, 2005; Zhang et al, 2007). The limitation of any pooling 
strategy is the potential to mask expressions in subsets of the 
general population; however, the subsequent targets will be more 
generalisable to the broader population. To prevent spurious 
samples affecting this analysis, we carried out a stringent inclusion 
criteria based on RNA quality before inclusion in our study and 
only included the 12 samples that met our criteria, out of the 50 we 
screened. Furthermore, to ensure subsequent targets were repre- 
sentative, we carried out validation in unpooled selected samples 
and found the results were representative of individual samples and 
not a result of outliers within the pool. Finally, we were able to 
demonstrate the significance of these targets in external data 
cohorts by demonstrating that these same genes were important 
within previous transcriptome studies of SCs {n = 7 benign and 
n = S malignant) and clinical prostate cancer (n = 4 benign, n = 5 
primary prostate cancer and n = 4 metastatic cancers) (Varambally 
et al, 2005; Birnie et al, 2008). 

With this target discovery approach, we demonstrate that 
specific SC pathways were still maintained in early and advanced 
prostate cancer. Specifically, the identification of a novel deficiency 
in testosterone catabolism within normal human SCs, in addition 
to CSCs and advanced prostate cancer, highlights mechanisms that 
support a persistent 'addiction' of prostate cancer cells to AR 
ligands in CRPC (Reid et al, 2010; Williamson et al, 2012). 
Although previous studies of testosterone glucuronidation have 
concentrated on the role of class II UGT enzymes on androgen 
responsiveness within prostate cancer, it is both class I and class II 
UGT enzymes that are required to fully di-glucuronidate 
testosterone for cellular export (Murai et al, 2006). In particular, 
UGT1A8 is the only UDP glucuronidation enzyme uniquely 
capable of di-glucuronidation, which is required to fully catabolise 
DHT (Murai et al, 2006) and has a specific activity against C19 
steroids (androgens), that is, 2- to 20-fold greater than most class II 
UGT enzymes (Tukey and Strassburg, 2000). 

We have previously demonstrated that normal SCs of the 
prostate express low levels of active AR, but that differentiated cells 
express higher levels of the receptor with higher AR activity. In the 
normal setting, despite a deficiency in UGTIA enzymes, it appears 
the receptor is tightly regulated in SCs (Williamson et al, 2012). 
It is in the cancer setting, where AR function is already overactive 
and being actively driven (Feldman and Feldman, 2001; Attard 
et al, 2009; Bonkhoff and Berges, 2010), that deficiency in UGTIA 
seems to contribute to pathogenesis, with low expression of the 
enzymes correlating with poor prognosis and downregulation 
in vitro seen to contribute to AR activation. Deficiency in UGTIA 
enzymes provides an additional mechanism to maintain AR 
signalling in CRPC, contributing to recognised phenomena of 
androgen hypersensitivity, AR amplification and local testosterone 
production (Feldman and Feldman, 2001). Attempts were made to 



determine whether these cells had increased testosterone levels 
using mass spectroscopy; however, because of in vitro culture 
conditions with large quantities of steroids that can passively cross 
the cell surface, the system was saturated while preparing 
specimens, and we were unable to determine intracellular levels 
with accurate sensitivity. Nevertheless, PSA expression remains an 
excellent proxy for AR pathway upregulation in our data. 

Studies of UGTAl enzymes have been limited to date, in part 
because of the similarities of each enzyme to each other. Future 
studies characterising mechanisms of UGTAl enzyme gene 
silencing would be of therapeutic interest in prostate cancer. Our 
work supports models of CRPC where there is a persistent reUance 
of prostate cancer cells on AR signalling and AR ligands, 
reinforcing a mechanistic role for newer therapies that target 
testosterone and highlights a novel androgen catabolism pathway 
for target development. 
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